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1. Introduction 



Newtonian infinite volume grand canonical Gibbs measures on Riemannian man- 
ifolds have been studied for more than a decade [TJ [21 [31 H] ■ Recently, in the 
context of general relativity, the finite volume canonical ensemble was developed 
and compared to Newtonian analogs for an ideal gas of test particles, provided 
that the container of gas is within the range of validity of Fermi coordinates for 
an observer whose worldline follows a timelike Killing vector, and the container 
is stationary relative to the observer [5]. 

Here we consider grand canonical ensembles, with the assumption that the test 
particles constituting the gas or fluid do not affect the background metric ten- 
sor. One might expect that the theory of specifications and Gibbs states for 
Riemannian manifolds developed in [TJ [21 [3] , and the references therein, could 
be applied in a straightforward manner to Lorentzian space-time manifolds, thus 
providing a framework for general relativistic statistical mechanics. However, 
there are several impediments. First, space-time and energy-momentum coor- 
dinates are not easily decoupled as in the Newtonian case. One must work from 
the start with the cotangent bundle of space-time. In some spacetimes there 
is a nonzero energy of the vacuum that must be accounted for in a thermody- 
namic formalism. There are also ambiguities associated with the notion of a 
thermodynamic limit, since space-time can be foliated into space slices along an 
observer's timelike path in many different ways. In addition, the notion of po- 
tential energy, which is a feature of central importance in statistical mechanics, 
normally has no place in general relativity. These issues are considered in the 
sequel. 

For the purposes of statistical mechanics, our requirement that the spacetime 
possess a timelike Killing field is natural. The Killing vector field, tangent to 
the timelike path of an observer, determines a time coordinate for the system, 
and thus an energy component of the four-momentum. As a consequence, com- 
parisons to corresponding non relativistic statistical mechanical formulas are 
possible in some space-times for certain observers. Moreover, this choice of 
four-velocity forces the particle system to evolve in such a way that the geome- 
try of spacetime is unchanging along its worldsurface (since the Lie derivative of 
the metric vanishes in that direction), making it plausible on physical grounds 
for the system to reach equilibrium. 

The significance of the frame of reference in relativistic statistical mechanics was 
discussed in [5]. Even in Minkowski space-time, there is no Lorentz invariant 
thermal state for an ideal gas system. A thermal state is in equilibrium only in 
a preferred Lorentz frame in which the container of gas is at rest, and therefore 
breaks Lorentz invariance [5] . The natural generalization for statistical mechan- 
ics in a space-time with a timelike Killing field is an orthonormal tetrad in which 
the Killing vector determines the time axis. 
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In Sect. 2 of this paper, we develop a non rotating, orthonormal "Fermi- Walker- 
Killing" coordinate system in which the Killing vector serves as the time axis. 
The measure on phase space based on these coordinates satisfies a Liouville the- 
orem and is invariant under space coordinate transformations. Sect. 3 applies 
the grand canonical formalism already available for Riemannian manifolds to 
the Fermi surface defined in the previous section, and establishes notation. Sect. 
4 provides a general relativistic ideal gas law. The term "ideal gas" is some- 
what misleading in the context of general relativity. This is because a volume 
of gas subject to no forces is still affected by the curvature of space-time, and 
this corresponds to a Newtonian gas subject to gravitational, tidal, and in some 
instances "centrifugal forces," but otherwise "ideal." 

In Sect. 5, we define "Newtonian limit" as the limit of statistical mechanical 
quantities as the speed of light c — > oo for a class of spacetimes whose timelike 
Killing fields have a certain dependence on c. We prove convergence of relativis- 
tic Gibbs states evaluated on cylinder sets to their Newtonian counterparts for 
the ideal gas case, and define associated Newtonian pressures. This Newtonian 
limit is illustrated in detail in Sect. 6 for a gas of particles following a circular 
orbit in Kerr space-time. In Sect. 7, we prove uniqueness of the infinite volume 
Gibbs state in anti-de Sitter space for a class of equilibrium interaction poten- 
tials and calculate the infinite volume relativistic pressure for an ideal gas of 
test particles. Concluding remarks are given in Sect. 8. 

2. Particle systems on space-time manifolds 

Let (M,g) be a smooth four-dimensional Lorentzian manifold. The Levi-Civita 
connection is denoted by V, and throughout we use the sign conventions of 
Misner, Thorne and Wheeler [TJ. Let cr(r) be a a timelike path parameterized 
by proper time r with unit tangent vector eo(r). The path, <t(t), represents the 
worldline of an observer, who (in principle) may make measurements. 

Measurements in a gravitational field are most easily interpreted through the 
use of a system of locally inertial coordinates. For the observer cr(r), Fermi- 
Walker coordinates provide such a system. A Fermi- Walker coordinate frame is 
nonrotating in the sense of Newtonian mechanics and is realized physically as 
a system of gyroscopes [TJ [HI [SI [ID]- Expressed in these coordinates, the met- 
ric along the path is Minkowskian, with first order corrections away from the 
path that depend only on the four-acceleration of the observer. If the path is 
geodesic, the coordinates are commonly referred to as Fermi or Fermi normal 
coordinates, and the metric is Minkowskian to first order near the path with 
second order corrections due only to curvature of the space-time [7j . A partial 
listing of the numerous applications of Fermi- Walker coordinates is given in [IT] . 
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A Fermi- Walker coordinate system for a may be constructed as follows. A 
vector field v in M. is said to be Fermi- Walker transported along a if v satisfies 
the Fermi- Walker equations given in coordinate form by, 

F^) = v^iny = o. (i) 

where Q a p = a a up — u a ap, and a a is the four-acceleration along a. A Fermi- 
Walker coordinate system along a is determined by an orthonormal tetrad of 
vectors, eo(r), e2(r), es{r) Fermi- Walker transported along a. Fermi- 

Walker coordinates 3 relative to this tetrad are defined by, 

x° ( ex P CT (r)(^e 3 (r)) = r 

(2) 

x k (exp CT(r) (A J e J (r)) = A fc , 

where here and below, Greek indices run over 0,1,2,3 and Latin over 1,2,3. 
The exponential map, exp p (i7), denotes the evaluation at affine parameter 1 of 
the geodesic starting at the point p in the space-time, with initial derivative v, 
and it is assumed that the \ 3 are sufficiently small so that the exponential maps 
in Eq.([2]) arc defined. From the theory of differential equations, a solution to 
the geodesic equations depends smoothly on its initial data so it follows from 
Eq.® that Fermi- Walker coordinates are smooth. Moreover, it follows from 
[T2"] that there exists a neighborhood of a on which the map \ — (x°, x 1 , x 2 , x 3 ) 
is a diffeomorphism onto an open set in R 4 and hence a coordinate chart. Gen- 
eral formulas in the form of Taylor expansions for coordinate transformations 
to and from Fermi- Walker coordinates, are given in |13j . 

Assume now that if is a timelike Killing vector field in a neighborhood of a and 
the tangent vector to a is K, i.e., eo = K. The vector field K is the infinitesimal 
generator of a local one-parameter group <j) s of diffeomorphisms. The function 
<fi s is the flow with tangent vector K at each point, and <t(t) = <fi T (0, 0, 0, 0) 
(where (0,0,0,0) is the origin in Fermi- Walker coordinates). 

Following [5], define a diffeomorphism x _1 from a sufficiently small neighbor- 
hood of the origin in M 4 to a neighborhood U of <r(0) in Ai by 

x-V,* 1 ^ 2 ,* 3 ) = toOcHo,* 1 ,* 2 ,*?)) (3) 

Then x = (x° , x 1 , x 2 , x 3 ) is a coordinate system on U which, as in [5], we refer 
to as Fermi- Walker-Killing coordinates. It follows that, 
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(0,x\x 2 ,i 3 ) = (0,z\x 2 ,x 3 ) 



_d_ 

dx 1 



- A 
x°=o dx l 
d _^ 



x°=0 



(4) 



Thus, in Fermi- Walker-Killing coordinates, the time coordinate x° is the pa- 
rameter of the flow generated by the Killing field K with initial positions of the 
form (0 3 ) in Fermi- Walker coordinates. 

In the case of Fermi- Walker coordinates we designate momentum form coordi- 
nates as p a and in the case of Fermi- Walker-Killing coordinates, the momentum 
form coordinates will be designated as p a . Thus, coordinates of the cotangent 
bundle of M. will be represented as {x a ,pp} or {x a ,pp}. The metric components 
in Fermi- Walker coordinates are designated as g a p and in Fermi- Walker-Killing 
coordinates as g a p- Note that g a/ 3 does not depend on x° because of Eqs.((4]). 

The state of a single particle consists of its four space-time coordinates together 
with its four-momentum coordinates {x Q ,p^}, but the one-particle Hamiltonian 
satisfies, 

H = \g ap P a P^~^f, (5) 

where, for timelike geodesies, m is the rest mass of the particle, and c is the speed 
of light. Eq.© just states the standard fact from relativity that the square of 
the "norm" of the four- velocity of a timelike particle is — c 2 . An arbitrary choice 
of pi,p 2 ,p3, determines p by, 



_ -g^Pi + V(9° 3 9 0k ~ 9 00 g' k )PjPk - 5 00 m 2 c 2 

Po - ^ , (6) 

so that Eq.© holds. 
Let ip : M -> K by, 

tp(x) = g(exp~{ Q) x, cr'(0)). (7) 

The space slice, X, at time coordinate r = = x°, orthogonal to the observer's 
four- velocity (along u(t)) is given by, 

X = ^- x (0). (8) 
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The restriction 3 g of the space-time metric g to X makes (X, 3 g) a Riemannian 
manifold, sometimes referred to as a Fermi surface or Landau submanifold, c.f., 
for example, |14j and references therein. 

Remark 1. The Fermi surface X defines the collection of space-time points 
that are simultaneous with c(0) in a natural way relative to the observer a. 
Hyper surf aces parallel to X foliate a (possibly global) neighborhood of a , and 
define a notion of simultaneity. This foliation has been used to study geometri- 
cally defined relative velocities of distant objects and time dependent diameters 
of Robertson- Walker cosmologies U51 \16j/ 

Thus, using Eq.©, the phase space for a single particle at time coordinate 
a; = is determined by the Fermi surface, X, together with the associated 
momentum coordinates, i.e., the cotangent bundle of the Fermi surface at fixed 
time coordinate x° . Let P be the sub bundle of the cotangent bundle of space- 
time with each three dimensional fiber determined by ©• The proof of the 
following proposition is given in |17j . 

Proposition 1. The volume 1-form G) on P given by, 

Q = -^-dx A dx 1 A dx 2 A dx 3 A dpi A dp 2 A dp 3 , (9) 
is invariant under all coordinate transformations. 

As discussed in [5], phase space for a single particle is determined by the Fermi 
surface, at fixed time coordinate, orthogonal to the observer's four-velocity 
(along ct(t)), along with the associated momentum coordinates, i.e., the cotan- 
gent bundle of the Fermi surface at fixed time coordinate x°. The appropriate 
volume form is given by the interior product i(md / dr)u) of the four-momentum 
vector md/dr = p a d/dx a with Cj. Then, 

mi(d/dT)u} = dx 1 A dx 2 A dx 3 A dpx A dp 2 A dp 3 + (dx° A tp), (10) 

where tp is a five-form. Since dx° = on vectors on phase space (with fixed 
time coordinate x°), the restriction of mi(d/dT)Q to the one-particle (six- 
dimensional) phase space is, 

mi(d/dT)ui = dx 1 A dx 2 A dx 3 A dpi A dp 2 A dp 3 . (11) 

It follows from Proposition[T]that the 6-form given by Eq. (|ll[) is invariant under 
coordinate changes of the space variables with x° fixed. Thus, 
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dxdp = dx 1 A dx 2 A dx 3 A dpi A dp2 A <ip3 
= dx 1 A dx 2 A cfe 3 A dp% A dp 2 A dp 3 
= dx 1 dx 2 dx 3 dp\dp2dpz 
= dxdp 

Physically this means that the calculations that follow are independent of the 
choice of the orthonormal triad ei(0), e2(0), es(0). 

In what follows, we consider a particle system in contact with a heat bath, 
and consequently in thermal equilibrium. Following the usual convention, let 
f3 = 1/kT, where T is the temperature of the gas in some volume A, and k is 
Boltzmann's constant. 



Remark 2. An alternative convention for temperature may be used in what 
follows, consistent with an idea of Tolman's in the context of relativistic ther- 
modynamics flSf . At any location within the gas, one may identify an observer 
with four-velocity K/\\K\\. The energy of a particle with momentum p mea- 
sured by that observer is —K a p a /\\K\\. Then, defining a position dependent 
equilibrium temperature T by T = T/\\K\\ = T / ^/—goo coincides with Tolman's 
formula and results in expressions equivalent to what we find in the sequel. 

It was shown in [5] that the energy of a single particle may be naturally defined 
as —K a p a — mc 2 , and that for a suitable volume A (see below) the one particle 
partition function is given by, 



Q l(/3 ,A)= * / f e^^^p = ^f [ ^x, 
{2-kHY J a J r3 {2TTh) A J A 7 

(13) 

where -^2(7) is the modified Bessel function of the second kind evaluated at 7, 
and, 

j = j(x\x 2 ,x 3 ) ee \\K\\f3mc. (14) 
Remark 3. It is easy to verify that, 

£d/dso (e^ Rap " +mc2 Ump) = C K (e^ Rap « +mc2 U5tdp) = 0, (15) 

where C denotes Lie derivative. Analogous to the development of nonrelativistic 
statistical mechanics, this invariance together with Eq. (|12[) is a partial justifica- 
tion for our choice of phase space measure. 
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We take the phase space volume form for n particles to be the product of n 
copies of Eq.(fT2j). one copy for each particle, which we denote as, 



dx n dp n = dx\dx\dx\dpndpi2dpi3 . . . dx n dx n dx n dp n idp n 2dp n 3 



(16) 



It was shown in [5] that this measure satisfies a Liouville theorem. We now as- 
sume that the n-particle equilibrium distribution is determined by the following 
partition function: 



?nC8,A,«) 



47r(mc)V mc 

(27Tft) 3 



e 0C£?=l K a p la -V A (^\s)) dy .n dv n 

71 / _ \ 

^VA(x"|si I TT ^(7i 



(17) 



A" 



n 



dx E 



7i 



where s denotes a configuration of particles outside of A that influences the 
equilibrium configuration of particles x n in A (s can be the empty configura- 
tion), Va(x"|s) = Va(x{, x\ , x\ , . . . x\, x„, x^ \s) is a many body potential energy 
function whose form and restrictions are given in the following section, and 
7i = j(xl,xf,xf). We note that this assumption excludes interaction functions 
of both position and momentum coordinates, a significant restriction. However, 
a similar assumption for the case of special relativity was made in [TI5] . We em- 
phasize that we are not asserting that potential energy functions determine (via 
a Hamiltonian) the dynamics of general relativistic particle systems; instead, we 
assume only that equilibrium behavior is governed by expressions of the form 
of Eq.([l7|. 



3. Grand canonical ensemble on the Fermi space slice 

Because the integrals in Eq. pTl) factor as a product of configuration integrals 
and momentum integrals (which can be evaluated explicity), we focus on the 
former. To proceed further, it is convenient to work with the grand canonical 
ensemble. A grand canonical formalism for particle systems on Riemannian 
manifolds is described in [TJ [3J [3] and references therein, briefly summarized 
here in a way that is useful for our purposes. 

Let B(X) denote the c-algebra of Borel sets on X and B C (X) represent the 
collection of sets in B{X) with compact closure. Let T x be the set of particle 
configurations on X , i.e., 

T x = {x C X : \x n K\ < oo for any compact K c X}, (18) 
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where \A\ denotes cardinality of the set A. For A G B C (X), define AT A : — > 
No, the set of natural numbers (including zero) by, 

N A (x) = \xD A|. (19) 

Let B(Tx) be the c-algebra generated by all such functions. For A C X, let 
ia = iI1A and, 

r A = {x e r x ■. x x/A = 0}. (20) 

For n = 0, 1,2, let, 

r A n) ,{ier A : |z| = n} . (21) 
There is a natural bijection, 

A n /S n ->r£°, (22) 
where 5* n is the permutation group over {1, n} and, 

A n = {(x u ...,x n )eA n :x i ^x j ifi^j}. (23) 

If A G B C (X) is open, this correspondence determines a locally compact, metriz- 
able topology on and then, 

r A = Q if > (24) 

is equipped with the topology of disjoint union and hence a Borel er-algebra, 
S(Fa). Then (Tx,B(Tx)) is the projective limit of the measurable spaces 
(Fa,S(Fa)) as A increases to X, pQ. 

For A G B C {X), let 6 A = cr{jV A , :A'cA, A' G B C (X)}. The cr-algebras B(T A ) 
and are u-isomorphic. 

Definition 1. Let A G B C (X). A set A G B A C B(Fx) is called a cylinder set. 

Remark 4. If A is a cylinder set, let A A = {x A : x G A}. Then 1a(x) — 
1a a {xa)- To see this, letC = {A c Tx ■ ^a(x) = 1a(^aUsx/a) f° r a ^ s € ^x}- 
Then 1a{x) — 1.a a (x a ) holds if and only if A G C. Let T> be the class of the 
sets of the form {x : N Al (x) — n%,N\ 2 (x) = U2, ■-, N Ak {x) = n^}, for some 
nonnegative integers k, tlx, 71%, rifc, and Ai, A2, A& C A. Then T> is closed 
under intersection and T> C C. Moreover, C has the following properties (i) 
T x G C, (ii) if A,B G C and A C B then B/A G C and (Hi) if A n G C 
and A n f A, then A G C. By Dynkin's Monotone Class Theorem, C contains 
o-(V)=B A . 
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We next define a Poisson measure on X, based on the measure dx, determined 
by Eg .IP]). Let A e B C (X) be open and let si : A" ->■ r£° by s£((xi, x„)) = 
{xi, ...,x n }. Then dx™ o (s^) _1 is a measure on r£° which we again denote by 
fix" (distinguished from Eq. (fTB| by context). For the case n = 0, e?x°(0) = 1. 
For nonnegative activity z, and A € 23 C (X), we can define (un normalized) 
Poisson measure on T\ by, 

oo n 

v A (dx) = V ^-dx". (25) 



n! 

n=0 



For each A € B C (X), let Va be a S(Fa) measurable interaction potential of the 
form, 



n=0 yCx,\y\=n 

for finite configuration x, where <p n : ->1U {+oo}, inf„^ cj)(y) > — oo. For 
the case n — 0, we define 4>q{0) = Va{0) = |A|p vac , where |A| is the volume of 
A determined by the Riemannian metric 3 g on X and p vac is the energy density 
of the vacuum given by 

Pvac = — ^, (27) 

07TG 

where A is the cosmological constant (possibly zero) and G is Newton's universal 
gravitational constant. 



Definition 2. An interaction potential of the form given by Eq. \2d>\ has finite 
range if for each n > 2, there exists R > such that <f> n (x\, x n ) = whenever 
the proper distance ( as determined by the Riemannian metric 3 g ) from Xi to Xj 
exceeds R for some i and j . 

Though not essential, we assume henceforth that each (f> n has finite range for 
n ^ 3, so that only the two-body interactions might have infinite range. In 
order to develop notational consistency with Eq. (|17|) , we also define a tempera- 
ture dependent, one-body potential (or equivalently, a temperature independent 
external field) by, 

-/3<MzW 3 ) - ^2(7^ W 3 )) /0fi ^ 
e — — -: — g — r\ — • 

J(X L ,X Z , X 6 ) 

For configurations x C A and s C X define, 
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VA(x\s) = \A\p vac + Y2 J2 Mv)> (29) 

n=l yCxll(snX/A) 
\y\=n 

if the right hand side converges absolutely, and +00 otherwise. For A £ B C (X) 
and a boundary configuration s, define the partition function by, 

Z A (s) = Z A (s,f3,z)= [ e-P v ^v A {dx), (30) 
Jr A 

where we assume that a factor 47r(mc) 3 e /3mc2 /(27r/i) 3 is included in the activity 
z through Eq.(|2"5"|) so that, 

where the parameter /1 is chemical potential. With this identification, Eg. (13111) 
may be rewritten as, 



Z A (s,(3,z) =J2ni Qn(PA,s), (32) 

n=0 

where Qo(j3,A, s) = exp(— (3\A\p vac ) and for n > 1, 

Q n ((3,A,s) = / e-^WV. (33) 

Following the grand canonical statistical mechanical prescription, the connection 
to thermodynamics is given by the following definitions. 

Definition 3. The finite volume pressure P A (/3, z) for A € B C {X) is given by, 

RP ( o v \ log^A(0) . . 

PPa[P,z) = , (34) 

where the volume \A\ determined by the Riemannian metric 3 g on X. 

Definition 4. For the case that X is unbounded, we define the infinite volume 
pressure P = P(/3, z) by, 

W ,,) = hmi^M, (35) 

n^oo |A n | 

provided the limit exists. Here A n is the ball of radius n centered at the origin 
of Fermi- Walker- Killing coordinates. 
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Following [3J, Eqs. ([23|) - (|30p determine a local specification II. For any s € , 

A e S C (X), and A € S(X), let, 



Ha(A s) = lz M<°°W [ 1 A ( X u s x/A )e-W*\>W(dx). (36) 
We define a Gibbs state via the DLR equations [3J . 

Definition 5. Assume that X is unbounded. A probability measure fx on 
(Tx , B(Tx)) is a Gibbs state if 

[ n A (A,s)fi(ds) = li{A) 

JVx 

for all A e B(T X ) and A e B C (X). 

Remark 5. For the case of pair potentials, conditions have been given for the 
existence of Gibbs states in J3]/; see also 0]/. Note also that the specification given 
by Eq. (I36p , and hence any Gibbs state, is invariant under changes of nonzero 
energy density of the vacuum 4>q (see ). However, the pressure is not. 

4. Ideal Gas Laws 

In the context of general relativity, we define an ideal gas to be the ensemble 
determined by a potential V satisfying <f> n = for all n > 2. Note, however, 
that (j>i (see Eq. (|28|) ). and (f>o, which is proportional to the vacuum energy 
given by Eq. (|27|) . is not necessarily zero. Using the definitions of the preceding 
sections, a derivation of a finite volume general relativistic ideal gas law for 
an observer following a path tangent to a timelike Killing vector, using the 
grandcanoncial ensemble is straightforward. By Def. [3J 

P^^ + ^^M*. (37, 
The expected number of particles, (N A ), in A is given by, 

(38) 

Combining Eqs. (J3TJ) and (|38]) then yields, 

P\ = -Pvac + "^p fcT = ^Pvac + P\kT, (39) 

where p A is the number density of particles in A. In the case that p vac = 0, this 
is the standard form of the ideal gas law, with the volume | A| determined by the 
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induced metric 3 g on the Fermi surface for the observer. The result is well-known 
for the case of Minkowski space-time. As in the classical (Newtonian) case, it 
is easily shown (using Eq. (j3"5)l ) that the particle number is Poisson distributed. 
Thus, if PrA(iV) represents the probability that N particles are in A, then, 

Pr A (i\0 = e-<"A>J^. (40) 

It follows from Eqs. (l3"51) . (13T)1) . and (|4U| that the thermodynamic behavior of 
an ideal gas, and in particular the particle density, is determined by the met- 
ric of the space-time through the norm of the timelike Killing field, (see 
Eq . (fl4"|) ) . and the volume |A| (determined by 3 g). Thus in principle, the geom- 
etry of space-time may be studied through the behavior of an ideal gas within it. 

By way of illustration, we compare the thermodynamic behavior of an ideal gas 
in de Sitter space and the Einstein static universe. It was shown in that the 
metric for the Einstein static universes in polar form using Fermi coordinates 
for any geodesic observer (with path a(t) = (t, 0,0,0) in these coordinates) is 
given by, 

Einstein static: ds 2 = -c 2 dt 2 + dp 2 + R 2 sin 2 (d9 2 + sin 2 6 d<t> 2 ), (41) 

where p < nR and R = 1/VX is the radius of the universe. It was also shown in 
[11] that the metric for de Sitter space-time, in polar form of Fermi coordinates 
for any geodesic observer is given by, 



de Sitter: ds 2 = -c 2 cos 2 (ap) dt 2 + dp 2 + sm2 fr) U6 2 + sin 2 6 d<j> 2 ), (42) 

a 

where a = yfX/Z and p < ir/2a so that the coordinates cover the space-time up 
to the cosmological horizon of the observer. For both space-times, K — d/dt is 
a timelike Killing vector field tangent to the respective timelike geodesies of the 
observers. Then from Eq. (fl4")) , we find for these space-times, 

Einstein static: ^(x 1 ,x 2 , x 3 ) — (3mc 2 , (43) 
as distinguished from, 

de Sitter: "/(x 1 , x 2 , x 3 ) = f3mc 2 cos (ap) , (44) 



where in terms of Cartesian Fermi coordinates, p — y/ (x 1 ) 2 + (a; 2 ) 2 + (a; 3 ) 2 . 
The expected number of particles for each space-time is determined by Eq.(|38j). 
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and respective volumes are calculated from Eqs.pTj) and P^j) after setting dt = 
0. For volumes determined by a given set of Fermi coordinates for the two 
space-times, Eg. (159)) and (1271) yield, 



PA=P= + Pmass— , (45) 

where p mass = mp\ is the mass density of the ideal gas when m is the mass 
of each gas particle, and these densities for the two space-times depend on the 
respective values of 7, and will in general be different. So the relativistic ther- 
modynamic behavior of an ideal gas distinguishes the two space-times for any 
common value of A. 



Remark 6. It is interesting to note that the equation of state for the perfect 
fluid of the Einstein static universe is P — 0, with zero temperature, i.e., the 
fluid is "dust, " and the mass density p determined by the Einstein field equations 



<46) 

If we set Pmass for the ideal gas of test particles equal to p of the perfect fluid 
and require P = 0, the temperature is given by, 

kT=^f, (47) 

and approaches zero only as m — > 0. Thus, the thermodynamic behavior of the 
ideal gas models that of the perfect fluid (i.e. dust) only in the limiting case that 
mass of each test particle of the gas goes to zero. 



5. Newtonian limits 

For a timelike path a tangent to a Killing vector if in a given space-time, we 
have defined through Eqs. ([57j) - P0|) the relativistic grand canonical pressure 
of an ideal gas. Under appropriate conditions, there is associated with such 
a statistical mechanical system an analogous Newtonian (i.e., nonrelativistic) 
expression for the pressure of the gas. It should be noted that the term "ideal 
gas" is somewhat misleading in the context of general relativity. This is be- 
cause a volume of gas subject to no forces is still affected by the curvature of 
space-time, and this corresponds to a Newtonian gas subject to gravitational, 
tidal, and in some instances "centrifugal forces," depending on the path of the 
observer. This associated Newtonian partition function and finite volume Gibbs 
state for an ideal gas is obtained as a limit as c — > 00. 
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Definition 6. For a timelike Killing field K tangent to the path a parameterized 
by proper time r, define a dimensionless function of the space coordinates on 
the Fermi surface, a = a(x , x 2 , x 3 ), by \\K\\ = ac. 

Throughout this section we let, y = 1/c so that the Newtonian limit is obtained 
as y — > + , and we assume that a can be extended to a smooth function of 
y,x l ,x 2 ,x 3 including at y — (here and below we suppress the dependence 
of a on x l ,x 2 ,x 3 ). We refer to this function as a(y) and further assume that 
a(0) = 1 and a'(0) = (see the following section). Then, 



cm* 



M y 3 



-a"(0)y 2 + O(y 3 



(48) 



where < yo < y. Note that o/'(0) is a function of x 1 , x 2 , x 3 . Below we will iden- 
tify |ma"(0) as the Newtonian, nonrelativistic potential energy U(x 1 , x 2 , x 3 ) of 
a test particle of mass m with coordinates (0 ) due to the gravitational 

field. This potential energy function U{x l ,x 2 ,x 3 ) is normalized (by an addi- 
tive constant) so that f/(0, 0, 0) = 0, i.e., the potential energy on the timelike 
path a is zero. This is because K is the four velocity on a so a(y) = \\K\\/c is 
identically 1 as a function of y = 1/c on cr, forcing a"(0) = there. 

Definition 7. Let A € B C (X). For a particle system following a timelike path 
a, we define the following Newtonian statistical mechanical analogs. 



(a) The Newtonian activity paramater, z^ ewt is given by, 



z Newt 



2irh 2 /3 

where [i is the chemical potential ( as in Eq 



(49) 



(b) Zfj ewt {N) denotes the partition function for a non relativistic ideal gas 
subject to an external field with potential energy function ^ma"(0) and 
activity z^ewt so that, 



n\ J \ 



= exp 



e -^ma"(0) dx 



z New t I e~^ ma "^ rfx 

A 



(50) 



(c) Let Hj[ ewt (A, s) denote the specification associated with the one-body New- 
tonian potential ^f3ma" (0) and partition function Z^ ewt ( A) given bu Ea. §3G\i . 
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Theorem 1. Let a be a timelike path with a satisfying i?o. (|48[) . Assume that 
A e B C (X). Then 

(a) If Z\ is the partition function for an ideal gas for an observer in a space- 
time with p vac = (i.e. with zero cosmological constant), then, 



lim Z A = Z Newt (A), 

c—>oo 

where as above c is the speed of light, 
(b) If A iE B a then, 

lim n A (A,s) = Il% ewt (A, S ) 



(51) 



(52) 



( where both sides of the equat 



ion are in 



dependent of s). 



Proof, a) The asymptotic behavior of -^2(7) for large argument is given by (see 

Ml 



-7 



(53) 



7 y 2 7 

Since K is timelike, ||2?|| is nonzero and bounded below on A. It follows that 
7 = ||A"||/3toc — > 00 uniformly on A as c —¥ 00. By ()53|) . given any e > 0, there 
exists S > such that 1/c = y < S implies, 



1 



< e. 



(54) 



Multiplying by z, substituting Eqs. (f3lj) . (|49|) . and 7 = a(3mc 2 , and rearranging 
terms yields, 



-/30i 



,/3mc ! (l-a) 



z Nowt" 



v3/2 



/3mc (1— a) 



,3/2 



(55) 



Since a is smooth and the closure of A is compact, it follows from Eq. (|4"5)l that 
a — > 1 and (3mc 2 (l — a) — > — f3ma" (0)/2, both uniformly on A as y — > 0. It 
then follows from Eq. ([55l) and the triangle inequality that, 



ze-P^ -> z Ncwt e-5/3-«"(o) ; 



(56) 



uniformly on A as y — > 0. Integrating both sides of Eq. (|5"6")) then gives, 

lim zQ x = z Ncwt f e -^ ma "^ rfx. (57) 



16 



Thus, 



lim Z\ = lim exp[zQi] = exp 



ZNewt / e 2 
A 



i/3ma"(0) dx 



Z Newt (A). (58) 



b) Let A e Sa be a cylinder set, and let A a = {x\ : x £ A}. By Remark [H 

l A (a;) = Wa*). (59) 

Since A A C T A we may write A A = U%L A% where = A^nT^. From Eqs. 
©and®, 

n A (u, s ) = — £-f/ iA S (^,...,^)n( ze " Wl(a:<) ) dxn - ( 6 °) 

A n=0 n ' t=l 

By Ea.(|56p. there exists Co independent of (a; 1 , cc 2 , x 3 ) such that 

ze~^ < z N cwte-^ mQ " (0) + 1, (61) 

for all c > cq and all (j; 1 , x 2 , a; 3 ) g A. Since a"(0) is bounded on A, there is a 
constant C such that 

oo n 

£ - / U X (a;i,..., a; „)n(^Newte-^ mQ "W + l)dx n 



n=0 

■DC 



(62) 



< Y^-. / U X (x 1 ,...,x„)C"dx"<e c ^ ldx <oo. 

Therefore part (b) follows from the Dominated Convergence Theorem and Eq. ([58| 

□ 

Corollary 1. Assume that the Fermi surface X is unbounded and let /j. be the 
unique Gibbs state for the ideal gas. Let fi^ewt oe the unique Newtonian Gibbs 
state for the specification {Hj[ ewt }. Then for any cylinder set A, 

lim fi{A) = n Newt (A). (63) 

Proof. Existence and uniqueness of the Gibbs state for an ideal gas follows from 
Kolmogorov's theorem for projective limit spaces (see e.g. [3] p. 21 or [21] Chap 
V, Theorem 3.2). Since A is a cylinder set, there exists A £ B C (X) such that 
Ha(A,s) is independent of s. It follows from Eq.© and Theorem [T] that, 



lim (i(A) = lim U A (A,s)= lim U^ ewt {A, s) = fts^A). (64) 



c— >oo 

□ 
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Definition 8. Let Pa((3,z), given by Eg. (|37p . be the pressure of an ideal gas 
according to a timelike observer satisfying the hypotheses of Theorem [TJa. We 
define the associated Newtonian pressure PNewt(P, z, A) to be, 

(3P Newt ((3, z, A) = Qf 1 dx^) log ^^(A). (65) 
In i/ie case that p vac ^ 0, we define the associated Newtonian pressure by, 

/3P Newt (f3,z,A) = -p vac + (^J ldxj \ogZ Newt (A). (66) 

We note that in the preceding definition, the role of volume A is restricted solely 
to the identification of (constant) limits of integration for the volume integral. 
In the relativistic context, these limits of integration are determined by proper 
lengths of the dimensions of the container of gas, and in the Newtonian limit 
the limits of integration are absolute length measurements. 

Our inclusion of the vacuum energy density term in Eq. (|66[) is motivated in 
part by generalized Buchdahl inequalities that establish a minimum density of 
matter in space-times with a positive cosmological constant [35]. Evidently, a 
collection of particles falling below a certain critical density will be pushed apart 
by a positive vacuum energy. 



6. Kerr space-time 

In this section we calculate a(y) and the limiting Newtonian potential energy 
function, U(x 1 , x 2 7 x 3 ) = ma"(0)/2 for the case of circular geodesic orbit in 
the equatorial plane of Kerr space-time. The Kerr metric in Boyer-Lindquist 
coordinates is given by (recall that y = 1/c), 



where, 



, 2 / 2y 2 GMr\ 1 , 2 Ay 2 GM ar sin 2 6 , S . 2/1Ij2 

ds 2 = - 1 \^dt 2 ^ dtdcj) + sin 2 Odcf 

\ p z ) y z p z p l 

2 

+ ^dr 2 + p 2 d9 2 , 



E= (r 2 +y 2 a 2 ) -j/VAf 



(67) 



p 2 = r 2 + y 2 a 2 cos 2 6, 

A = r 2 -2GMy 2 r + y 2 a 2 , (68) 
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and where G is the gravitational constant, M is mass, and a is the angular 
momentum per unit mass, and —GMy ^ a ^ GMy. 

Below we will need notation for A, and S, evaluated at the specific coordinates 
r = tq, and 9 = n/2. For that purpose, we define, 

A = r 2 -2GMy 2 r + y 2 a 2 , 
s o = (To + y a ) -y a A . 

The circular orbit in the equatorial plane (see, e.g., [23], [24]) with d(j)/dt > 
at radial coordinate rp, is given by, 



tr(t)= U,r ,7r/2, — , (70) 

When a > 0, this orbit is co-rotational and when a < the orbit is retrograde. 
All circular orbits in the equatorial plane are stable for tq > 9GMy 2 , indepen- 
dent of a, but the co-rotational circular orbit is stable for smaller values, even 
inside the ergosphere at r = 2GMy 2 when a is sufficiently close to GMy 2 (cf. 
[23]). 

The Killing field tangent to the path (where it is the four velocity) is given by, 

g= ^VGM^ M vmy (71) 

where, 



V = y2y 2 aVGMr { i/2> +4-3y 2 GMrl (72) 

Using Eqs. flBTl) and ([7T]) . a(y) = yy / —K a K a may be found in Boyer-Lindquist 
coordinates. Straightforward calculations show that a(0) = 1 and a'(0) = 
and, 



1 •» . . GMm GMmr 2 sin 2 9 3GM m , . 
-ma"(fi) = + — . (73) 

2 w r 2rl 2r 

This expression may be interpreted as Newtonian potential energy of a single 
particle in a container of particles in circular orbit with angular velocity, 
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where ro is the radius of the orbit. The first term on the right hand side of 
Eq. (|73p is the gravitational potential energy, and the third term is an additive 
constant that forces a"(0) = when r = ro and 6 = ir/2, i.e., at the origin of 
coordinates for the rotating container. The second term is the centrifugal poten- 
tial energy, more readily recognized when expressed in cylindrical coordinates. 
To that end, let p = r sin#, with <f> the azimuthal angle, and let z measure linear 
distance along the axis of rotation. The magnitude of angular momentum of 
a uniformly rotating particle of mass m with cylindrical coordinates (p, <f>, z) is 
I = mp 2 4>, and the centrifugal potential energy is then given by, 

I 2 GMmp 2 _ GMmr 2 sin 2 6 

The minus signs in Eq. (|75[) take into account the direction of force, away from 
the central mass. 



In order to compute a"(0) in Fermi coordinates, we select the following tetrad 
vectors in the tangent space at cr(0). 



eo 



' y 2 ay/GM + rf /2) VGm \ 

v ' ' ' v r 



ei = (°'^r' '°) ' 

e 2 = (o,0,i-,oJ, (76) 

(y 2 VGM{y 2 a 2 + r 2 -2y 2 aVGM^) y 2 VGM(a - 2^) + r^ 3/2) \ 
e 3 — ~ n- ,U,U,- 



This tetrad may be extended via parallel transport to the entire circular orbit 
given by Eq. ([70)) . but we need these tetrad vectors only at cr(0). 

Fermi coordinates relative to these coordinate axes (i.e., the above tetrad) may 
be calculated using Eq.(27) of [13]. The result for the Boyer-Lindquist coordi- 
nates r and 9 to second order expressed in the Fermi space coordinates 
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at (a; = t = t = 0) is, 



7T X 2 y/~KoX 1 X 2 

2 r r2 



+ 



(77) 



'•o 



/A^x 1 , (r y 2 GM - y 2 a 2 ) (x 1 



TO 



2r 3 



Ao^ 2 ) 2 (r -y 2 GM)(x 3 Y 



2r 3 



2r 2 



(78) 



Now, calculating a(y) = '—K a K a in Boycr-Lindquist coordinates, substitut- 
ing for r and 9 using Eqs. ([771) and ([781 gives, 



. . 1 GMy 2 (r 2 +3a 2 y 2 ~ Aay 2 ^GM) (x 2 f 

a (y) = 1 + kzt- 



2r 2 V 2 

3 GMy 2 A (x 1 " 



2 °~" (79) 



2r 2 V 2 



0(3) 



Computing the second derivative with respect to y at y = yields, 



1 „, , / GMm GMmx 1 GMm (2 (i 1 ) 2 - (i 2 ) 2 _ (r 3 ) 2 ) 
_„,„»(„) . _____ + __ V __ 



GMm GMmx 1 GMmUx 1 ) + (x 3 ) . i 3GM 



+ + Ho L + + 0(3) 



1 -m 



2r ' r 2 2rg J ' 2r 

(80) 

Eq. (|80l) may be compared term- by-term with Eq. (l73[) . The expression in the 
first pair of parentheses on the right hand side of Eq. ((80|) is the Taylor expan- 
sion to second order of the gravitational potential, which is the first term on 
the right hand side of Eq. (j?3"|) . The second terms in both equations are related 
analogously. At the point a(0) in the orbit, the Cartesian variable x 1 in Eq. flBTfl) 
measures (Newtonian) distance from the origin of coordinates in the radial di- 
rection away from the central mass, x 2 measures distance in the "z direction" 
parallel to the axis of rotation, and x 3 measures distance in the tangential di- 
rection, parallel to the motion of the container of gas in orbit. (However, these 
orientations do not hold at other parts of the orbit since the coordinate axes are 
nonrotating.) Eq. (|80[) may obviously be simplified to yield the potential energy 
function, 
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U(x\x 2 ,x 3 ) = -ma"(0) = -±-L + -4--L + (3). (81) 

z zr zr 



Combining Eq. ([8T|) with Eqs.(|50j) and (|65|) immediately yields the Newtonian 
pressure (with y = 1/c) for a gas in circular orbit around a central mass to 
which the relativistic counterpart given by Eq. ([37j) may be compared. 



7. Anti-de Sitter space 

In this section, we prove uniqueness of the infinite volume Gibbs state in anti-de 
Sitter space for a class of equilibrium interaction potentials and calculate the 
infinite volume relativistic pressure for an ideal gas of test particles. 

Explicit transformation formulas to and from Fermi coordinates {x°, x , x 2 , x 3 } 
for a geodesic observer in (the covering space for) anti-de Sitter space-time are 
given in The geodesic path of the Fermi observer in these coordinates is 
<r(t) = (t, 0,0,0). Fermi coordinates for this space-time are global, i.e., one co- 
ordinate patch covers the entire space-time. 

Under the change of coordinates, x° = t, x 1 — ps'm9 cos0, x 2 = p sin 9 sin <j>, 
x 3 = pcos9, the metric becomes diagonal. In these "polar-Fermi coordinates," 
for any appropriate range of the angular coordinates, the line element becomes, 

ds 2 = -c 2 cosh 2 (ap) dt 2 + dp 2 + Slnh2 9 (ap) (d9 2 + sin 2 6d0 2 ), (82) 

where a = -y/|A|/3 and A < is the cosmological constant. For this space-time, 
K = d/dt is a timelike Killing vector field tangent to a so that eo = K. From 
Eqs.CEl and (|82]l. 

-fix 1 , x 2 , x 3 ) = \\K\\/3mc = I3mc 2 cosh(ap), (83) 

where p = y^x 1 ) 2 + (x 2 ) 2 + (a; 3 ) 2 . The Fermi surface, X, given by Eq.(0), is 
hyperbolic space. 

Because a single chart, with domain R 3 , covers the entire Fermi surface for this 
space-time, each configuration of particles on the Fermi surface corresponds to 
a configuration on R 3 , and Gibbs states on the Fermi surface are in a natural 
one-to-one correspondence with Gibbs states of configurations on R 3 with Pois- 
son measure of Eq. (|25|) based on Lebesgue measure dx. For an elaboration, see 
[U [21 E] - In this context, and noting Remark 03 we may regard V as an interac- 
tion on configurations on R 3 and apply the following special case of Theorem 1 
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of [2S] to the present circumstances. 



Theorem 2. Suppose there exists an increasing sequence of bounded Borel sets 
{A;.} whose union is R 3 such that: 

(1) (j> n (x) = 0, for any n > 2 and x — (xi, x n ) such that ifl / ^ 
x n (K 3 /A fc+ i) for some k; 

(%) J2k I su Ps z A k (s, P, z)\~ diverges, where A k = A k+1 /A k . 
Then there is at most one Gibbs state for V, p,z. 

Assume that V is a superstable interaction and has the form given by Eqs. ([28)) 
and (|2T)1) . Assume further that V has finite range and satisfies the condition, 

oo 

V A (x 1 ,...,x n \s) = Y^ <l>n{y)> -Bn, (84) 

n=2 yCxU(snX/A) 
\y\=n 

for some B > 0, all boundary configurations s, and all configurations (xi, x n ). 
Observe that V is the same as V except that V does not include the vacuum 
energy 4>o or the one-body potential function, (f>\. The potential V satisfies 
Eq. (|54l if it is positive or satisfies a hard-core condition (see [25]). 



Theorem 3. Let V be a potential energy function with finite range R satisfying 
i?g. ([84[) on the Fermi surface X in anti-de Sitter space. Then there is at most 
one Gibbs state on (Tx,B(Tx)) f or V-,P-, Z f or a ^ choices of P and z. 

Proof. Let Afe C X be the sphere of radius kR centered at the origin of coordi- 
nates in the chart for X, where R is as in Def. [2] Then condition (1) of Theorem 
[2]is satisfied by {A^}. To see that condition (2) of Theorem[5]is also satisfied, 
let A k be defined as in Theorem [5] and observe by Eg. (154]) that, 



sup Z Ak (s,P, z) < / exp 

* JTa, 



PBn-pJ2M^) 



y Ak {dx). (85) 



It follows from Egs.([53"]) and ([53]) that when p = ^/(x 1 ) 2 + (x 2 ) 2 + (x 3 ) 2 is 
sufficiently large, 



P4>l (x 1 ,x 2 ,x 3 ) ^ /3mc 2 cosh(ap) 



(86) 



Then for large k, 
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sup Z A Js, P, z) < exp\3Bn - /3ncosh(kR)]—r \A k \ n , (87) 
a ^ n\ 

n— 

where \A k \ = J Ak dx 1 dx 2 dx 3 , i.e., the Lebesgue measure of the image of A k in 
the Fermi coordinate chart. Hence, |j4fe| < C(kR) 2 for a constant C. Thus, 

supZa»(«,jM < exp[zC(i?/fc) 2 e /3s -^ cosh(fcfl) ], (88) 

8 

and condition (2) of Theorem [2] is clearly satisfied. □ 

Corollary 2. The pressure of an ideal gas of test particles in anti-de Sitter 
space is given by, 

'~i3? < 89 > 

where A < is the cosmological constant for anti-de Sitter space. 

Proof. With the same notation as in the proof of Theorem |5J using Eq. (157)1 and 
changing to spherical coordinates, we find, 

47rz f k K 2 (f3mc 2 cosh(ap)) 
f3mc 2 \A k \ J cosh(ap) 

The result now follows from Eqs.flSSJ), dUJ, and flS5J) 



P Ak 08, z) = - Pmc + I V H " P 2 dp. (90) 



□ 

Remark 7. The infinite volume pressure of an ideal gas of test particles in anti- 
de Sitter space is independent of temperature and chemical activity, and depends 
only on the magnitude of the cosmological constant. There is no corresponding 
Newtonian analog, as developed in Sect. 5, because j?q.(|48p does not hold for 
the Killing vector d/dt. 

8. Concluding Remarks 

For space-times with timelike Killing fields, we have developed a grand canonical 
formalism for statistical mechanical systems of test particles. The thermody- 
namic behavior of an ideal gas was shown in Sect. 4 to be determined by the 
metric of the space-time through the norm of the timelike Killing field, and the 
volume of the particle system as determined by projection 3 g of the metric on 
the Fermi space slice of the timelike observer. 

We proved uniqueness of the Gibbs state for a class of interaction potentials in 
anti-de Sitter space, and found the infinite volume pressure for an ideal gas of 
test particles. For the case of an ideal gas in space-times where Eq.(|48[) is sat- 
isfied, we introduced a notion of Newtonian limit of Gibbs states and of grand 
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canonical pressures. The physical legitimacy of the grand canonical formalism 
was demonstrated in Sect. 6, where it was shown that the Newtonian limit of a 
relativistic particle system following a circular geodesic orbit in Kerr space-time 
results in exactly what is predicted by classical mechanics. 

Directions for further research include possible generalizations of the statistical 
mechanical formalism which assume only the existence of approximate timelike 
Killing vector fields or conformal timelike Killing fields, relaxing our assumption 
of the existence of a timelike Killing vector field. In addition, alternative notions 
of simultaneity associated with coordinate systems other than Fermi- Walker- 
Killing coordinates might be considered, such as standard curvature-normalized 
coordinates on Robertson- Walker cosmologies for observers co-moving with the 
Hubble flow, or optical coordinates (in which events are "simultaneous" if they 
are visible to the observer a at the same proper time coordinate of a). In this 
way, the effect of the curvature of space on statistical mechanical behavior of 
test particles might be analyzed for realistic circumstances. 

Acknowledgment. The authors wish to thank Peter Collas for helpful discus- 
sions. 
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